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ABSTRACT 

We report the discovery of two transiting extrasolar planets by the HATSouth survey. HATS-9b 
orbits an old (10.8 d= 1.5 Gyr) V=13.3 G dwarf star, with a period P & 1.9153 d. The host star has a 
mass of 1.03 M©, radius of 1.503R© and effective temperature 5366 =b 70 K. The planetary companion 
has a mass of 0.837Mj, and radius of 1.065 Rj yielding a mean density of 0.85gem -3 . HATS-lOb 
orbits a V=13.1 G dwarf star, with a period P ~ 3.3128 d. The host star has a mass of 1.1 M©, radius 
of 1.11 Rq and effective temperature 5880 =b 120 K. The planetary companion has a mass of 0.53 Mj, 
and radius of 0.97 Rj yielding a mean density of 0.7gem -3 . Both planets are compact in comparison 
with planets receiving similar irradiation from their host stars, and lie in the nominal coordinates of 
Field 7 of K2 but only HATS-9b falls on working silicon. Future characterization of HATS-9b with 
the exquisite photometric precision of the Kepler telescope may provide measurements of its reflected 
light signature. 

Subject headings: planetary systems — stars: individual ( HATS-9, GSC 6305-02502, HATS-10, 

GSC 6311-00085 ) techniques: spectroscopic, photometric 


1. INTRODUCTION 

Our current understanding of the structure and or¬ 
bital evolution of extrasolar giant planets has been, to 
a large degree, informed by the characterization of tran¬ 
siting planetary systems. Besides the determination of 
the planet radius, true mass, and bulk density, follow-up 
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studies of transiting extrasolar planets (TEPs) allow the 
extraction of valuable information, like the spin-orbit an¬ 
gle and the properties and composition of the planetary 
atmospheres, that are not be easily recovered unless the 
orbital plane is favorably oriented such that the planet 
eclipses its host star. 

Detections of giant TEPs, mostly driven b y transiting 


groun d based surve ys like SuperWASP (Pollacco et al., 
20061 and HATNet ( Bakos et al.||2004 ), have revealed a' 


large number of systems in the region of parameter space 
with R p >0.8Rj, M p >0.4Mj , P<5 d and FGK-type host 
stars. The measured properties of these systems, coupled 
with subsequent follow-up studies, have been fundamen¬ 
tal for testing formation and interior models of these gi¬ 
ant planets, which are known as hot Jupiters. 

New gro un d-base d transiting surveys like HATSouth 
(Bakos et al. 2013) have been designed with the goal 
ol expanding the parameter space of well characterized 
TEPs by detecting planets with smaller radii ( R p < 
OARj) and/or longer periods (P > 10 days). In the 
process of searching for these kinds of planets, new 
hot Jupiters are detected which contribute to enlarging 
the sample of known systems. Even though many hot 
Jupiters are already known, more are still needed to make 
headway into understanding their physical properties, 
e.g. a firm understanding of the mechanism that causes 
some hot Jupiter s to h ave infl ated radii (e.g. HAT-P-32b 
and HAT-P-33b, Hartman et al. 120111). 

New planet discoveries around bright stars, accessible 
by follow-up facilities are especially valuable given the 
wealth of detailed studies that they can be subject to. In¬ 
deed, some of the most analyzed and characterized giant 
TEPs are three planets (TrES-2b, HAT-P-7b and H AT- 


P-llb) that were detected by g r ound-based survey s (Pal 
. 12010 O’Donovan et al. 20061) 
AyA’s Kepler mission (Borucki 


et al. 2008 Bakos et al 
and later o bserved by N ASA’s Kepler mission (Borucki 
et al. 2010). Even though the primary goal oi Kepler 
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satellite was the detection of planets near the habitable 
zone for estimating their frequency and distribution in 
our galaxy, the high photometric precision of Kepler al¬ 
lowed very detailed studies of the small population of 
giant planets on close orbits around moderately bright 
stars (V < 14) that fell in its field of view. 

Kepler was able to detect secondary tra nsits and phase 
varia tions on TrES-2b and HAT-P-7b (Esteves et al. 
2013) which were useful in the study of their atmospher- 
ical properties, such as the determination of the geomet¬ 
ric albedos and planetary phase curve offsets. Doppler 
beaming and ellipsoidal variations measured with Kepler 
also constrained the mass of those planets. In the case 
of HAT-P-llb, Kepler observations were useful in char¬ 
acterizing the activity of the K-type host star; and the 
analysis of crossing stellar spots allowed the determine 
tion of the spin-orbit misalignment of this sy stem (Dem- 
ing et al.|2011; Sanchis-Ojeda & Winn|2011 j. Simultane¬ 
ous observations of the transits of HAT-P-llb by Kepler 
and Spitzer allowed also the detection of w ater vapor in 
the a tmosphere of this Neptune-size planet (Frame et al. 
2014|. Estimation of the planetary physical parameters 


depend strongly on the estimated stellar properties. In 
this regard, Kepler was also able to measure model inde¬ 
pendent stellar properties by the use o f asterose ismology 
on the thre e mentioned systems (Christensen-Dalsgaard 
etahpOi O). 


After the failure of two of its reaction wheels, the Ke¬ 
pler satellite is still working, but with a new observation 
strategy and a photometric precision within a factor of 
2 of the nominal Kepler mission performance (e.g., Van- 


derburg & Johnson 


grain et al. 


Mackey et al.|2015| | (Jrossfield et al. |2015 ). This 


2015 Foreman- 


new mis¬ 
sion concept, called K2 (Howell et ai.|2U14 ), will observe 
10 fields, each for a period oFapproximateiy 70 days, and 
some of these fields lie in the Southern hemisphere. One 
of the limitations of K2 is that the number of stars that 
can be monitored in each held is substantially lower than 
for the original Kepler mission. For this reason, the pre¬ 
selection of targets, based on ground-based observations 
of K2 fields is especially important for an efficient use of 
the satellite. 

In this work we present the discovery of HATS-9b and 
HATS-lOb, two hot Jupiters discovered by the HAT- 
South survey which are located in the nominal coordi¬ 
nates of Field 7 of K2 mission. In Section [2] we sum¬ 
marize the observations that allowed the discovery and 
confirmation of these planets. In Section [3] we show the 
global analysis of the spectroscopic and photometric data 
that confirmed the planetary nature of the transiting can¬ 
didates and also rejected blend scenarios that can mimic 
the photometric and radial velocity signals. Our findings 
are discussed in Section |4] 

2. OBSERVATIONS 
2.1. Photometric detection 

HATS-9 and HATS-10 were identified as transiting 
planetary host candidates after obtaining ^10000 im¬ 
ages of the same held with three stations on the three 
HATSouth observing sites. The number of photomet¬ 
ric observations that were taken for each star on each 
of the HATSouth stations is indicated in Table [l] where 
it can be seen that in both cases ~45% of the observa¬ 


tions came from the HATSouth station located at Las 
Campanas Observatory (LCO). 

The HATSouth observations consist of four-minute 
Sloan r-band exposures obtained with 24 Takahashi E180 
astrographs (18cm aperture) coupled to Apogee 4Kx4K 
U16M ALTA CCDs. Readout times are of the order of 
one minute which results in a cadence of about 5 min¬ 
utes. Detailed descriptions of the image processing steps 
and the candidate id entification procedu res o f HATSouth 
data c an be found in Bakos et al. (2013) and Penev et al. 
(20131. Briefly, after applying aperture photometry on 
the images, the light curves generated are detrended us¬ 
ing external parameter decorrelation (E PD) and trend 
filtering algorithm (TFA Kovacs et al. 2005). Periodic 
transits on the detrended light curves are then s earched 
usin g Box-f itted Least Squares (BLS) algorithm (Kovacs 
et~aL|[2002 


Figure [l] shows the phase-folded detection light curves 
of HATS-9b and HATS-lOb, where a clear ~10 mrnag 
flat-bottom transit can be observed in both cases. 

2.2. Spectroscopic Observations 

Transit-like light curves can be produced by different 
configurations of stellar binaries. Spectroscopic observa¬ 
tions are required to reject false positives and to obtain 
the orbital parameters and masses of the true planets. 
Due to the great number of HATSouth candidates and 
the limited available observing time on spectroscopic fa¬ 
cilities, this follow-up is performed in a two-step proce¬ 
dure as we now describe. All spectroscopic observations 
are summarized in Table [2j 

First, initial spectra are acquired (with either low res¬ 
olution, or low S/N) to make a rough estimation of 
the stellar parameters, identifying spectra composed of 
more than one star, and measuring RV variations pro¬ 
duced by stellar mass companions. HATS-9 was ob¬ 
served with WIFeS (Dopita et al. 2007) on the ANU 2.3m 
telescope, obtaining 'Jog*=5821 ± 300 K, logg*=3.9 ± 
0.3 [Fe/H]=0.5 ± 0.5; and with ARCES on the APO 
3.5m obtaining T eff *=5692 ± 50 K, logg*=4.14 ± 0.1 
[Fe/H]=0.5 ± 0.08. Both estimates of stellar parameters 
were consistent with a G-type dwarf, but the sub-solar 
surface gravity value points towards a slightly evolved 
system. Details on the observing strategy, reduction 
methods and t he processing of th e spectra for WIFeS 
can be found in Bayliss et al. ( 2013]) . The ARCES obser¬ 
vation was carried out using the 176 x 3'.'2 slit yielding an 
echelle spectrum with 107 orders covering the wavelength 
range 3200-10000A at a resolution of AA/A ~ 31500. A 
single ThAr lamp spectrum was obtained immediately 
following the science exposure with the telescope still 
pointed toward HATS-9. The science observation was 
reduced to a wavelength-calibrated spectrum using the 
standard IRAF echelle package^ and analyzed using_the 
Spectral Pa rameter Classification (SPC) program (Buch- 
have et al. 2012) to determine the radial velocity and 
stellar "atmospheric parameters. 

Reconnaissance spectroscopy was performed for 
HATS-10 using the echelle spectrograph mounted on 

1 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Phase-folded unbinned HATSouth light curves for HATS-9 (left) and HATS-10 (right). In each case we show two panels. The 
top panel shows the full light curve, while the bottom panel shows the light curve zoomed-in on the transit. The solid lines show the model 
fits to the light curves. The dark filled circles in the bottom panels show the light curves binned in phase with a bin size of 0.002. 


Table 1 

Summary of photometric observations 


Instrument/Field 3, 

Date(s) 

# Images 

Cadence b 

(sec) 

Filter 

Precision 0 

(mmag) 

HATS-9 

HS-1/G579 

2010 Mar-2011 Aug 

4317 

300 

r band 

6.9 

HS-3/G579 

2010 Mar-2011 Aug 

2138 

303 

r band 

7.6 

HS-5/G579 

2010 Sep-2011 Aug 

2784 

303 

r band 

6.9 

FTS 

2013 Apr 11 

134 

80 

i band 

1.4 

PEST 

2013 May 31 

186 

130 

Rc band 

3.4 

HATS-10 

HS-1/G579 

2009 Sep-2011 Aug 

4389 

301 

r band 

7.3 

HS-3/G579 

2010 Mar-2011 Aug 

2596 

303 

r band 

7.2 

HS-5/G579 

2011 Mar-2011 Aug 

3297 

303 

r band 

7.8 

CTIO 0.9m 

2012 Aug 29 

69 

213 

z band 

2.3 

FTS 

2013 Apr 05 

142 

63 

i band 

4.3 

GROND 

2013 Jun 14 

92 

156 

g band 

0.8 

GROND 

2013 Jun 14 

88 

156 

r band 

1.3 

GROND 

2013 Jun 14 

94 

156 

i band 

0.7 

GROND 

2013 Jun 14 

89 

156 

z band 

0.8 

PEST 

2013 Jun 27 

145 

130 

Rc band 

4.6 


a For HATSouth data we list the HATSouth unit and field name from which the observations are taken. HS-1 and -2 are located at Las Campanas 
Observatory in Chile, HS-3 and -4 are located at the H.E.S.S. site in Namibia, and HS-5 and -6 are located at Siding Spring Observatory in Australia. 
Each field corresponds to one of 838 fixed pointings used to cover the full A-n celestial sphere. All data from a given HATSouth field are reduced 
together, while detrending through External Parameter Decorrelation (EPD) is done independently for each unique field+unit combination. 
b The median time between consecutive images rounded to the nearest second. Due to weather, the day-night cycle, guiding and focus corrections, 
and other factors, the cadence is only approximately uniform over short timescales. 
c The RMS of the residuals from the best-fit model. 


the du Pont 2.5m telescope at Las Campanas Observa¬ 
tory. One observation using the 1" x 4" slit (AA/A ~ 
40000) was enough to confirm that HATS-10 has a sin¬ 
gle lined spectrum with the following stellar parameters: 
T e ff*=6100 ± 100 K, log=4.6 ± 0.5 [Fe/H]=0.0 ± 0.5, 
v sin *=5.0 ± 2.0 krn/s. This spectrum was reduced and 
analysed with an automated pipeline developed to deal 
with data coming from a host of different echelle spectro¬ 
graphs (Brahm et al. in prep.). The pipeline for DuPont 
is very similar to the ones we ha ve previou sly detailed 
for Coralie and FEROS data in Jordan et ah ( 2014). 

Once both candidates were identified as single-lined 
late-type dwarfs, spectra from high precision instruments 
were required to measure RV variations with high pre¬ 
cision (< 30 m/s) in order to measure the mass of the 
substellar companions and obtain the orbital parameters. 
HATS-9 and HATS-10 were observed several times with 


Coralie ( Queloz et al. 2001) on t he 1 2m Euler telescope, 
FEROS (Hauler & Pasquini 1998) on the 2.2m MPG tele¬ 
scope, and BPS on the 8m Subaru telescope (Noguchi] 


et al. [2002). Coralie and F EROS data were pr ocessed 
with the pipeline described in Jordan et al. (2014), where 
RV values are obtained using the cross correlation tech¬ 
nique against a binary mask and bisector span (BS) mea¬ 
surements a re com puted from the cross-correlation peak 
following Queloz et al. (2001). HDS RVs were measured 
using the procedure detailed in Sato et al. (2002, 2012) 
which are in turn based on the method of Butler et al. 
(1996) while BS values were obtained following Bakos] 


et al. (2007). 

Phased high-precision RV and BS measurements are 
shown for each system in Figure [2] and the data are 
listed in Tabic [3] Both candidates show RV variations in 
phase with photometric ephemeris, however for HATS-10 
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the residuals are higher than expected. This deviation 
can be partly explained by moonlight contamination in 
5 spectra acquired with Coralie in August 2013 which 
are marked with crosses in Figure [2l There are no sig¬ 
nificant correlations between RV and BS variations and 
thus we conclude the RV variations are not produced 
by stellar activity. The 95% confidence interval for the 
Pearson correlation coefficient between RV and BS was 
computed for both candidates using a bootstrap proce¬ 
dure. The confidence intervals are [-0.57,0.07] and [-0.43, 
0.37] for HATS-9 and HATS-10, respectively. The indi¬ 
vidual FEROS spectra were median combined for both 
candidates to perform a precise estimation of the stellar 
parameters. 


2.3. Photometric follow-up observations 

In order to confirm the occurrence of the transits and 
to better constrain the orbital and physical parameters 
of the companions, higher precision light curves for both 
candidates were acquired using several telescopes around 
the globe. Table [T| summarizes the key aspects of this 
photometric follow-up, including the dates of the obser¬ 
vations, the cadence and the filter. 

Two partial transits of HATS-9 were detected using 
the 0.3m Perth Exoplanet Survey Telescope (PEST) and 
the Spectral camera on the 2m Faulkes Telescope South 
(FTS), part of Las Cumbres Observatory Global Tele¬ 
scope (LCOGT). Results of these observations are pre¬ 
sented in Table 0] and shown in Figure [3] Two partial 
transits of HATsHO were observed with FTS and the 
CTIO 0.9m telescope. Another two full transits were 
measured with PEST and the GROND instrument on 
the MPG 2.2 m. These HATS-10 light curves are shown 
in FigureS] All the facilities used for high precision pho¬ 
tometric follow-up have been previously used by HAT- 
South; the instrument specifications, observation strate- 
gies and reduction p rocedures adopted c an be found in 
Bayliss et al.|(|2013|),|Zhou et al.|(|2014b|),|Hartman et al. 
120141) andiMohler-Fischer et al I d2013t) for FTS, PEST, 
CTIO 0.9m, and GRU1NU, respectively. 
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Table 2 

Summary of spectroscopy observations 


Instrument 

UT Date(s) 

# Spec. 

Res. 

AA/A/1000 

S/N Range 3 

7RV b 

(kms -1 ) 

RV Precision 0 
(ms- 1 * ) 

HATS-9 

APO 3.5 m/ARCES 

2012 Aug 25 

1 

31.5 

27 

-11.5 

500 

ANU 2.3 m/WiFeS 

2012 Sep 8 

1 

3 

140 



Euler 1.2 m/Coralie 

2012 Nov 6-10 

4 

60 

14-20 

-10.634 

37 

MPG 2.2 m/FEROS 

2012 Aug-2013 May 

9 

48 

32-76 

-10.653 

32 

Subaru 8m/HDS 

2012 Sep 19 

3 

60 

100-114 



Subaru 8m/HDS+l2 

2012 Sep 20-22 

9 

60 

60-100 


11 

HATS-10 

du Pont 2.5m/Echelle 

2013 Aug 21 

1 

40 

48 

-29.2 

500 

Euler 1.2 m/Coralie 

2012 Aug-2013 Aug 

12 

60 

17-23 

-28.131 

68 

MPG 2.2 m/FEROS 

2013 Mar-Jul 

5 

48 

29-85 

-28.044 

50 

Subaru 8m/HDS 

2012 Sep 22 

3 

60 

74-94 



Subaru 8m/HDS+l2 

2012 Sep 19-21 

9 

60 

41-99 


14 


a S/N per resolution element near 5180 A. 

k For Coralie and FEROS this is the systemic RV from fitting an orbit to the observations in Section l3.3l For ARCES and the du Pont Echelle 
it is the measured RV of the single observation. We do not provide this quantity for instruments for which only relative RVs are measured, or for 
WiFeS which was only used to measure stellar atmospheric parameters. 

c For High-precision RV observations included in the orbit determination this is the RV residuals from the best-fit orbit, for other instruments 
used for reconnaissance spectroscopy this is an estimate of the precision. 
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Figure 2. Phased high-precision RV measurements for HATS-9 (left), and HATS-10 (right) from HDS (filled circles), FEROS (open 
triangles), and Coralie (filled triangles). In each case we show three panels. The top panel shows the phased measurements together with 
our best-fit model (see Table [d]) for each system. Zero-phase corresponds to the time of mid-transit. The center-of-mass velocity has 
been subtracted. The second panel shows the velocity O — C residuals from the best fit. The error bars include the jitter terms listed 
in Table |6| added in quadrature to the formal errors for each instrument. The third panel shows the bisector spans (BS), with the mean 
value subtracted. Note the different vertical scales of the panels. RV measurements highly contaminated with moonlight are marked with 
crosses. 
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Table 3 

Relative radial velocities and bisector spans for HATS-9 and HATS-10. 


BJD 

RV a 

°RV b 

BS 

OBS 

Phase 

Instrument 

(2,456,000+) 

(ms- 1 ) 

(ms” 1 ) 

(ms -1 ) 

(ms" 1 ) 



HATS-9 

169.62456 

82.24 

28.00 

-126.0 

13.0 

0.686 

FEROS 

171.51911 

107.24 

26.00 

-66.0 

13.0 

0.675 

FEROS 

173.70707 

67.24 

37.00 

-93.0 

17.0 

0.817 

FEROS 

189.85616 



-3.6 

13.0 

0.249 

Subaru 

189.87089 



-2.3 

11.0 

0.257 

Subaru 

189.88561 



9.6 

9.2 

0.264 

Subaru 

190.84576 

129.27 

7.39 

6.2 

13.4 

0.766 

Subaru 

190.86048 

146.04 

8.17 

5.3 

10.0 

0.773 

Subaru 

190.87520 

134.32 

7.39 

-6.7 

13.6 

0.781 

Subaru 

191.84685 

-139.82 

13.91 

-9.6 

21.4 

0.288 

Subaru 

191.86157 

-133.96 

12.35 

6.0 

13.5 

0.296 

Subaru 

191.87633 

-112.04 

12.85 

19.1 

21.9 

0.304 

Subaru 

192.84272 

112.22 

8.57 

-3.7 

18.5 

0.808 

Subaru 

192.85745 

126.06 

10.89 

-28.3 

16.3 

0.816 

Subaru 

192.87217 

100.81 

9.05 

8.0 

8.8 

0.824 

Subaru 

205.55663 

-48.76 

22.00 

2.0 

11.0 

0.446 

FEROS 

213.50471 

74.24 

24.00 

6.0 

12.0 

0.596 

FEROS 

215.55235 

125.24 

30.00 

-31.0 

14.0 

0.665 

FEROS 

219.55921 

172.24 

23.00 

-3.0 

12.0 

0.757 

FEROS 

237.50625 

-67.84 

28.00 

41.0 

22.0 

0.128 

Coralie 

238.50417 

86.16 

31.00 

157.0 

24.0 

0.649 

Coralie 

239.50502 

-141.84 

36.00 

16.0 

24.0 

0.171 

Coralie 

241.53678 

-82.84 

64.00 

349.0 

32.0 

0.232 

Coralie 

424.89567 

64.24 

36.00 

66.0 

17.0 

0.965 

FEROS 

427.91875 

7.24 

45.00 

-20.0 

20.0 

0.544 

FEROS 

HATS-10 

160.60805 

-132.00 

28.00 

-68.0 

21.0 

0.267 

Coralie 

161.58511 

-77.00 

34.00 

29.0 

22.0 

0.561 

Coralie 

164.61796 

-98.00 

33.00 

12.0 

22.0 

0.477 

Coralie 

189.90597 

-27.97 

14.29 

-6.1 

27.3 

0.110 

Subaru 

189.92069 

-42.32 

22.31 



0.115 

Subaru 

189.93541 

-34.14 

35.70 



0.119 

Subaru 

190.89115 

-35.74 

13.59 

17.7 

38.3 

0.408 

Subaru 

190.90587 

-24.73 

14.10 

10.6 

18.6 

0.412 

Subaru 

190.92060 

-26.35 

14.32 

-17.0 

25.1 

0,416 

Subaru 

191.89225 

49.44 

23.96 

3.2 

15.3 

0.710 

Subaru 

191.90697 

62.06 

24.16 

0.2 

17.3 

0.714 

Subaru 

191.92170 

46.62 

38.06 



0.719 

Subaru 

192.88724 



-34.5 

31.0 

0.010 

Subaru 

192.89965 



-0.2 

23.1 

0.014 

Subaru 

192.91206 



25.9 

15.7 

0.018 

Subaru 

237.55535 

-46.00 

46.00 

123.0 

24.0 

0.493 

Coralie 

238.55388 

70.00 

37.00 

-24.0 

22.0 

0.795 

Coralie 

239.53128 

-108.00 

46.00 

-37.0 

24.0 

0.090 

Coralie 

241.51608 

177.00 

52.00 

67.0 

26.0 

0.689 

Coralie 

375.90634 

-44.76 

23.00 

16.0 

12.0 

0.255 

FEROS 

376.90643 

47.24 

24.00 

37.0 

12.0 

0.557 

FEROS 

377.90816 

109.24 

24.00 

56.0 

12.0 

0.860 

FEROS 

427.83173 

108.24 

51.00 

114.0 

22.0 

0.929 

FEROS 

491.79404 

-123.76 

20.00 

14.0 

10.0 

0.236 

FEROS 

524.51947 c 

-98.00 

29.00 

249.0 

21.0 

0.115 

Coralie 

524.59002° 

-59.00 

25.00 

196.0 

19.0 

0.136 

Coralie 

524.70382° 

-75.00 

29.00 

315.0 

21.0 

0.170 

Coralie 

525.53878° 

57.00 

39.00 

462.0 

24.0 

0.422 

Coralie 

525.65573° 

-103.00 

38.00 

459.0 

24.0 

0.458 

Coralie 


Note. — Note that for the iodine-free template exposures we do not measure the RV but do measure the BS and S index. Such template 
exposures can be distinguished by the missing RV value. The Subaru/HDS observations of HATS-10 without BS measurements have too low S/N 
in the l 2 -free blue spectral region to pass our quality threshold for calculating accurate BS values. 

a The zero-point of these velocities is arbitrary. An overall offset 7 re i fitted independently to the velocities from each instrument has been 
subtracted. 

k Internal errors excluding the component of astrophysical jitter considered in Section l3.3l 
k Coralie observations acquired in August 2013 were contaminated with moonlight. 
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Figure 3. Left: Unbinned transit light curves for HATS-9. The light curves have been corrected for quadratic trends in time fitted 
simultaneously with the transit model. The dates of the events, filters and instr ume nts used are indicated. The second curve is displaced 
vertically for clarity. Our best fit from the global modeling described in Section l3.3l is shown by the solid lines. Right: residuals from the 
fits are displayed in the same order as the left curves. The error bars represent the photon and background shot noise, plus the readout 
noise. 



Figure 4. Similar to Figure [3] here we show the follow-up light curves for HATS-10. 
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3. ANALYSIS 

3.1. Properties of the parent star 

We determine precise stellar parameters for HATS-9 
and HATS-10 using a new code called ZASPE (Zonal 
Atmospherical Stellar Parameter Estimator) on me¬ 
dian combined FEROS spectra. The detailed structure 
and performance of ZASPE will be presented elsewhere 
(Brahm et al. in preparation), but in summary ZASPE is 
a Python-based code that computes the y 2 between the 
observed spectr a and the PHOENIX grid of synthetic 
spectra (Husser et al. 2013) only on the most sensitive 
spectral zones to each stellar parameter. The optimal set 
of stellar parameters (T e ff*,fog<?*, [Fe/H] and usini) is 
found iteratively and the sensitive zones are determined 
in each iteration. One of the most novel features of ZA¬ 
SPE is that the errors on the stellar parameters are com¬ 
puted from the data itself and include the systematic 
mismatches between the observations and the best fitted 
model. We have validated the results of ZASPE against 
a set of stars w ith interferometrically determined stel¬ 
lar parameters (Boyajian et al. 2012) and which have 
publicly available EEROS spectra. Results of this com¬ 
parison are shown in Figure [5] The resulting parameters 
for HATS-9 are: T eff *=5363 ± 90 K, log 5 *=3.97 ± 0.2, 
[Fe/H]=0.33 ± 0.09, v sin i=4.67 ± 0.5 km/s; while for 
HATS-10 we get: T eff *=5974 ± 110 K. log=4.44 ± 
0.13, [Fe/H] =0.19 ± 0.07, nsini=5.66 ± 0.5 km/s. 

These sets of stellar parameters were refined using the 
information contained in the transit light-curves. The 
stellar mean density (p*) can be computed directly from 
one of the light-curve model parameters (a/R±) and the 
period and eccentricity of the orbit by using Kepler’s 
third law with only a slight dependence on the stel- 
l ar par a meters through the limb-darkening coefficients 
(Sozzetti et al. 12007). The spectroscopically determined 
'Jeff* and [Fe/HJ were coupl ed with p* an d Yonsei-Yale 
stellar evolution models (Y2;|Yi et al.|2001) to determine 
the stellar physical parameters (A*, iW* and the age of 
the star), which were used to compute a new and more 
precise estimation of logp* for HATS-9 (log p* =4.12 ± 
0.04) and HATS-10 (logp*=4.38 ± 0.03). A new set of 
T e ff*, [Fe/H], nsini was determined using ZASPE with 
logp* fixed to the precise values obtained by modeling 
the light curves, followed by a new estimation of p* and a 
new modeling of stellar isochrones. The new set of stellar 
parameters fixing logp*, which are the ones we adopted 
for further analysis, were consistent with the initial val¬ 
ues quoted in the previous paragraph and are listed in 
Table [5] where distances are determined by comparing 
the measured broad-band photometry listed in that ta¬ 
ble to the predicted magnitudes in each filter from the 
isochrones. We assu me a Ry = 3.1 extinction law from 
Cardelli et al.| (1989) to determine the extinction. The 
lcr and 2cr confidence ellipsoids in T e ff* and p* are plot¬ 
ted in Figure [6] for both planet hosts, along with the 
Y2 isochrones for the ZASPE determined [Fe/H]. We 
find that HATS-9 is a 1.030±0.039 M 0 , 1.503±g Rq, 
quite evolved (10.8 ± 1.5 Gyr) star, while HATS-10 is a 
1.101 ±0.054 M 0 , 1.105 +q oIo Rq main-sequence star. 

We attempted to measure the Lithium absorption line 
at 6707.8 A for testing the age estimation of HATS-9 but 
the quality of our spectra was only enough to rule out a 


strong absorption feature. 

3.2. Excluding blend scenarios 

In order to exclude blend scenarios we carried out a 
blend analys is of the observations following |Hartman| 


et al. (2012). For HATS-9 we find that scenarios in¬ 
volving blends between a stellar eclipsing binary and a 
foreground or background star can be ruled out with 
greater than 5<r confidence based on the photometric 
data alone. The primary constraint in this case is the 
lack of out-of-transit variations seen in the HATSouth 
light curve. Due to the short orbital period, the best- 
fit blend model which reproduces the shape of the tran¬ 
sit has a ~ 1 mrnag amplitude ellipsoidal variation, and 
a ~ 0.5mmag deep secondary eclipse, neither of which 
are detected in the HATSouth observations. Moreover 
the Subaru/HDS observations of HATS-9 show no sig¬ 
nificant bisector span variation (the RMS scatter of the 
BS measurements is 12ms -1 ) providing further evidence 
that the system is not a blended eclipsing binary. For 
HATS-10 the photometric observations can be fit by a 
G±M star eclipsing binary blended with another G star 
that is slightly brighter than the primary in the eclips¬ 
ing system. Based on the difference in y 2 , this model 
is indistinguishable from a single G star with a transit¬ 
ing planet. We simulated spectra for blend models that 
could plausibly fit the photometric observations, finding 
that it all cases the blended systems would have easily 
been detected as having composite spectra. They also 
would produce RV and BS variations of several kms -1 , 
whereas the observed RV variation is 67 ± 10ms _1 , and 
the Subaru/HDS BS scatter is only 18 ms -1 . We con¬ 
clude that neither HATS-9 nor HATS-10 is a blended 
eclipsing binary system. As is often the case, however, we 
are not able to rule out the possibility that either transit¬ 
ing planet system has a fainter stellar-mass companion. 
For both systems a stellar companion of any mass, up 
to the mass of the planet-hosting star, is possible. If a 
massive stellar companion is present in a given system, 
the true planet radius would be up to ~ 60% larger than 
inferred here. The planet mass would also be larger. 
High resolution adaptive optics imaging, and/or long¬ 
term RV observations are needed to determine whether 


either system has a stellar companion (e.g. 

Howell et al. 

2011 Horch et al. 2014 

Everett et al.||2015 



3.3. Global modeling of the data 
We modeled the HATSouth photometry, the follow- 


following Pal et al. 

(2008 

; 

Bakos et al. 

2010): 

Hart-| 

man et ah 

(2012 

). We tit M 

landel & Ago! 

(2002) 

tran- 


of the HATSouth transit depth as a result of blending 
from neighboring stars and over-correction by the trend¬ 
filtering method. For the follow-up light curves we in¬ 
clude a quadratic trend in time in our model for each 
event to correct for systematic errors in the photometry. 
We fit Keplerian orbits to the RV curves allowing the 
zero-point for each instrument to vary independently in 
the fit, and allowing for RV jitter which we we also vary 
as a free parameter for each instrument. 

We used a Differential Evolution Markov Chain Monte 


Carlo procedure (ter Braak 2006 Eastman et al. 2013) 
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Table 4 

Light curve data for HATS-9 and HATS-10. 


Object 8 - 

BJD b 

(2,400,000+) 

Mag c 

^Mag 

Mag(orig) d 

Filter 

Instrument 

HATS-9 

55744.07098 

-0.00037 

0.00552 

0.00000 

r 

HS 

HATS-9 

55749.81701 

-0.00018 

0.00572 

0.00000 

r 

HS 

HATS-9 

55780.46237 

0.00906 

0.00604 

0.00000 

r 

HS 

HATS-9 

55767.05534 

-0.01086 

0.00553 

0.00000 

r 

HS 

HATS-9 

55696.18926 

0.01787 

0.00581 

0.00000 

r 

HS 

HATS-9 

55657.88321 

-0.00168 

0.00549 

0.00000 

r 

HS 

HATS-9 

55726.83440 

0.00055 

0.00619 

0.00000 

r 

HS 

HATS-9 

55680.86732 

0.01840 

0.00534 

0.00000 

r 

HS 

HATS-9 

55788.12454 

0.00811 

0.00716 

0.00000 

r 

HS 

HATS-9 

55776.63287 

-0.00411 

0.00550 

0.00000 

r 

HS 


Note. — This table is available in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and 
content. 


a Either HATS-9, or HATS-10. 

k Barycentric Julian Date is computed directly from the UTC time without correction for leap seconds. 

c The out-of-transit level has been subtracted. For observations made with the HATSouth instruments (identified by “HS” in the “Instrument” 
column) these magnitudes have been corrected for trends using the EPD and TFA procedures applied prior to fitting the transit model. This 
procedure may lead to an artificial dilution in the transit depths. For HATS-9 our fit is consistent with no dilution, for HATS-10 the HATSouth 
transit depth is ~ 93% that of the true depth. For observations made with follow-up instruments (anything other than “HS” in the “Instrument” 
column), the magnitudes have been corrected for a quadratic trend in time fit simultaneously with the transit. 

d Raw magnitude values without correction for the quadratic trend in time. These are only reported for the follow-up observations. 




Figure 5. (Left): Comparison between T e ff values obtained with ZASPE and the ones derived from interferometric information. (Right): 
Comparison between log(g) values obtained with ZASPE and the ones derived from interferometric information. 
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Figure 6. Model isochrones from|Yi et al.|(|2001|) for the measured metallicities of HATS-9 (left) and HATS-10 (right). In each case we 
show models for ages of 0.2 Gyr and. l.U to l4.0Gyr in 1.0 Gyr increments (ages increasing from left to right). The adopted values of T e ff* 
and p* are shown together with their la and 2 a confidence ellipsoids. The initial values of T e ff* and p* from the first ZASPE and light 
curve analyzes are represented with a triangle. 


to explore the fitness landscape and to determine the posterior distribution of the parameters. 
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Table 5 

Stellar parameters for HATS-9 and HATS-10 



HATS-9 

HATS-10 


Parameter 

Value 

Value 

Source 

Astrometric properties and 

cross-identifications 



2MASS-ID. 

2MASS 19231442-2009587 

2MASS 19371363-2212161 


GSC-ID. 

GSC 6305-02502 

GSC 6311-00085 


R.A. (J2000). 

19 h 23 m i4 28s 

19 h 37 m 13.80s 

2MASS 

Dec. (J2000). 

—20°09'58.7" 

—22°12'16.1" 

2MASS 

MR.A. (masyr- 1 ) 

0.3 ±4.3 

3.1 ± 1.3 

UCAC4 

MDec. (masyr -1 ) 
Spectroscopic properties 

-1.9 ±2.8 

-3.2 ± 1.6 

UCAC4 

T eff * (K). 

5366 ± 70 

5880 ± 120 

ZASPE a 

[Fe/H]. 

0.340 ± 0.050 

0.15 ±0.10 

ZASPE 

v sin i (kms 1 ) .... 

4.58 ±0.90 

5.68 ±0.70 

ZASPE 

»mac (kms -1 ). 

4.6 

3.8 

Assumed b 

^mic (kms -1 ). 

1.0 

1.0 

Assumed 0 

7RV (kms -1 ). 

Photometric properties 

-10.644 ±0.013 

-28.088 ±0.024 

Coralie,FEROS 

B (mag). 

14.080 ±0.010 

13.820 ±0.010 

APASS d 

V (mag). 

13.276 ±0.010 

13.113 ±0.010 

APASS d 

9 (mag). 

13.629 ±0.010 

13.448 ±0.010 

APASS d 

r (mag). 

13.072 ±0.030 

12.967 ± 0.010 

APASS d 

* (mag). 

12.865 ±0.010 

12.781 ± 0.010 

APASS d 

J (mag). 

11.885 ±0.022 

11.866 ±0.024 

2MASS 

H (mag). 

11.558 ±0.027 

11.568 ±0.024 

2MASS 

K a (mag). 

Derived properties 

11.479 ±0.022 

11.511 ± 0.025 

2MASS 

M* (M 0 ). 

1.030 ±0.039 

1.101 ± 0.054 

YY±p*±ZASPE e 

R* (Re) . 

1 rr 0 o+0.101 

i.ouo_q 0 43 

1 inK+0-055 

1.1UO_ 0 040 

YY ±p* ±Z ASPE 

logfl* (cgs). 

4.095 ± 0.038 

4.392 ± 0.032 

YY ±p* ±ZASPE 

P * (gem -3 ). 

0 427+ 0 - 030 
u.'iz/ _o 070 

1 1 C)+0- 12 

YY±p* ±ZASPE f 

L* (£©). 

1 70+ 0 - 24 
i - ' u -0.16 

1.31 ± 0.18 

YY±p*±ZASPE 

M v (mag). 

4.33 ±0.15 

4.52 ±0.16 

YY±p*±ZASPE 

Mx (mag,ESO).... 

2.49 ±0.13 

3.05 ±0.10 

YY±p*±ZASPE 

Age (Gyr). 

10.8 ± 1.5 

3.3 ± 1.7 

YY±p*±ZASPE 

A v (mag). 

0.000 ±0.011 

0.112 ±0.075 

YY±p*± ZASPE 

Distance (pc). 

622t« 

496 ± 24 

YY±p*±ZASPE 


a ZASPE = **** routine for the analysis of high-resolution spectra (***cite***), applied to the FEROS spectra of HATS-9 and HATS-10. These 
parameters rely primarily on ZASPE, but have a small dependence also on the iterative analysis incorporating the isochrone search and global 

modeling of the data, as described in__ ___ 

b Computed following [Valenti &; Fischer) |2005|. 
c |Husser et al.| ( |2013| ). 

^ From APASS DR6 for HATS-9, HATS-10 as listed in the UCAC 4 catalog ( |Zacharias et al.|2012[ .. 
e YY+p*+ZASPE = Based on the YY isochrones ( |Yi et al.|2001 ), p* as a luminosity indicator, and the ZASPE results. 

^ In the case of p* the parameter is primarily determined from the global fit to the light curves and RV data. The value shown here also has a 
slight dependence on the stellar models and ZASPE parameters due to restricting the posterior distribution to combinations of p*+T e ff* + [Fe/H] 
that match to a YY stellar model. 


The resulting parameters for each system are listed in 
Table[6| HATS-9b has a radius of 1.065 ±0.098 i?j and a 
mass of 0.837 ± 0.029 Mj, while HATS-lOb has a radius 
of 0.9691^45 R i and a mass of 0.526 ± 0.081 M J- Both 
planets have bulk densities slightly lower than the one 
of Jupiter (0.85 ± 0.19 g cm -3 and 0.70 ± 0.15 g cm -3 , 
respectively) 

4. DISCUSSION 

We have presented the discovery of two new transiting 
planets which are shown on mass-radius and equilibrium 
temperature versus radius diagrams in FigureJT] From 
the mass-radius diagram, HATS-9b and HATSAOb can 
be classified as typical non-inflated hot Jupiters. HATS- 
9b is slightly less massive than Jupiter (0.84 Mj) and 
has almost the same radius. Its orbital period of P = 1.9 
days is rather short compared to the period distribution 


of known hot Jupiters. HATS-lOb has a mass in the 
range between Saturn and Jupiter (0.53 Mj), a radius 
consistent with that of Jupiter and a period of P = 3.3 
days, which is close to the mean period of known hot 
Jupiters. 

The equilibrium temperature versus radius diagram 
shows that both planets tend to depart from the known 
correlation between the planet radius and its deg ree of 
irra diatio n. This correlation, first proposed in |Guil-| 


lot (2005), indicates that the inflated radius of some 
hot Jupiters can be at least partially explained by the 
enhanced insolation from their parent star. HATS- 
9b has a moderately high equilibrium temperature 
(T eq =1823j(3gK) due to the small star-planet separation 
coupled to the large stellar radius, while HATS-lOb has 
a more typical equilibrium temperature for a hot Jupiter 
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Table 6 

Orbital and planetary parameters for HATS-9b and HATS-lOb 



HATS-9b 

HATS-lOb 

Parameter 

Value 

Value 

Light curve parameters 



P (days) . 

1.9153073 ±0.0000052 

3.3128460 ± 0.0000058 

T c (BJD) a . 

2456124.25896 ± 0.00086 

2456457.88193 ± 0.00022 

Ti 4 (days) a . 

0.1457 ± 0.0024 

0.1253 ±0.0011 

T 12 — T 34 (days) a . 

0.0106 ±0.0015 

0.01157 ± 0.00100 

a //?* . 

4 36+ 010 

O '70 + 0.29 
°* ,o -0.44 

C /R* b . 

14.84 ±0.26 

17.588 ± 0.067 

R P /R * . 

0.0725 ± 0.0041 

0.0903 ±0.0013 

Ir . 

0 071 + 0 - 112 

u - u ‘ 1 -0.050 

on 3+°° 87 

U.iid >_ 0 059 

b = a cos i / R+ . 

0.271°;“ 

0 34+0 - 11 

U.O^_o.io 

i (deg) . 

86.5±H 

87.79 ±0.72 

Limb-darkening coefficients c 
ci, g (linear term) . 


0.5380 

C 2 ,g (quadratic term) . 


0.2487 

ci,r . 

0.4688 

0.3459 

C2,r . 

0.2596 

0.3349 

ci ,i . 

0.3533 

0.2587 

C2,i . 

0.2892 

0.3388 

ci, z . 


0.1978 

C2,Z . 


0.3360 

Cl,i? . 

0.4369 

0.3216 

C 2 , R . 

0.2687 

0.3371 

RV parameters 

K (ms -1 ) . 

133.5 ±3.4 

67 ± 10 

e d . 

< 0.129 

< 0.501 

RV jitter HDS (ms" 1 ) e . 

0.1 ± 5.2 

0.00 ±0.53 

RV jitter FEROS (ms -1 ) .... 

0.0 ±1.7 

38 ± 28 

RV jitter Coralie (ms -1 ) . 

Planetary parameters 

0.0 ± 1.1 

45 ± 23 

M v (Mj) . 

0.837 ± 0.029 

0.526 ±0.081 

Rp (Rj) . 

1.065 ±0.098 

0 969 +0 061 
u-uoy _ 0 045 

C(M P , R p ) f . 

0.48 

0.02 

Pp (gem -3 ) . 

0.85 ±0.19 

0.70 ±0.15 

l°g 9p ( c g s ) . 

3.253 ±0.068 

3.140 ±0.082 

a (AU) . 

0.03048 ± 0.00038 

0.04491 ± 0.00074 

To q (K) . 

1 823±|| 

1407 ± 39 

08 . 

0.0460 ± 0.0039 

0.0440 ±0.0071 

logic (±> (cgs) h . 

9 397+ 0 049 
_ 0 033 

8.947 ± 0.047 


a Times are in Barycentric Julian Date calculated directly from UTC without correction for leap seconds. T c : Reference epoch of mid transit that 
minimizes the correlation with the orbital period. T 14 : total transit duration, time between first to last contact; T 12 = T 34 : ingress/egress time, 
time between first and second, or third and fourth contact. 

k Reciprocal of the half duration of the transit used as a jump paramete r ill. our MCMC analysis in place of a/ /?*. It is related to a/ R r by the 
expression C/R* — a/R* (2 tt(1 + e sin uj))/(P-\/l — b~\/ 1 — e 2 ) |Ebd<osem^h 010 1. 

c Values for a quadratic law, adopted from the tabulations by |^^^et[ “ fl200f [ according to the spectroscopic (ZASPE) parameters listed in Table [ 5 ] 
d As discussed in Section 13.31 the adopted parameters for all four systems are determined assuming circular orbits. We also list the 95% confidence 
upper limit on the eccentricity determined when cos oj and y^e since are allowed to vary in the fit. 

e Term added in quadrature to the formal RV uncertainties for each instrument. This is treated as a free parameter in the fitting routine. 

1 Correlation coefficient between the planetary mass M p and radius R p estimat ed from the posterior par ameter distribution. 
g The Safronov number is given by O — / ( 14.,;,,/Vi,,h ) 2 — ( 0 / Rp ) ( M s ,/M> ) (see [Hansen fc Barman|2007| ■ 

1 Incoming flux per unit surface area, averaged over the orbit. 


(T eq =1407 ± 39K). Accor ding; to the empirical relations 


proposed in Enoch et al. (2012), which give the radius 
of a giant planet from its equilibrium temperature and 
semi-major axis, HATS-9b and HATS-lOb should have 
radii of 1.36i?j and 1.22i?j, respectively. The observed 
radii are 3cr and 5<r below these values, which indicates 
that these planets are very compact given their irradi¬ 
ation levels and that thus additional variables must be 


10b may have significant amounts of heavy ele ments in 
their cores. According to the interior models of [Fortney 
et al. (2007), both planets will require a core mass of ~60 
to exp ain their radii based on their masses, stellar 
host masses and orbital periods for an age of 4.5 Gyr. 
This explanation can be further motivated by the rela¬ 
tively high metallicity of their parent stars (0.340 ±0.050 


responsible of setting the radii of short period giant plan- 

1 Guillot et al. 

2006 

Burrows et al. 2007 |Enoch et al. 

ets. 

2011 

2012) have proposed a correlation between the in- 


lerred core mass of giant planets and the metallicity of 
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the parent star. The principal idea behind the proposed 
correlation is that a more metal rich proto-planetary disk 
will be more efficient in creating massive cores following 
the core-accretion scenario of planetary formation. Even 
though this process is expected to occur in the forma¬ 
tion and migration steps, the final relation between the 
stellar metallicity and the radius of giant planets is not 
at all clear and other pheno mena can act in the oppo- 
site direction. As shown by Burrows et al. (2007), the 
presence of heavy elements in the atmosphere of young 
giant planets will increase its opacity, slowing the con¬ 
traction and making the planetary radius more inflated 
than expected. Moreover, the validity of the proposed 
correlation has been put into question by the analysis 


Zhou et al. (2014a) who find no significant correlation 
between R p and |Fe/H] for the complete sample of de¬ 
tected giant TEPs. 

The age of the system may be another important vari¬ 
able, since the radius of giant planets should undergo 
Kelvin-Helmholtz contraction as t hey age, control led by 
their upper radiative atmosphere (Hubbard 1977). Fig¬ 
ure [8] presents the mass-radius diagram of transiting 
hot Jupiters having similar insolation levels to HATS- 
9b (1750K < T e q < 1900K). This figure shows that in 
general the bloating of the atmosphere of strongly irradi¬ 
ated planets is prevented for more massive hot Jupiters. 
This correlation presents some outliers, with HATS-9b 
the most extreme one. A peculiarity of HATS-9b is the 
advanced age of the system (~11 Gyr) contrasted with 
the ages of the rest of the planets in Figure[8](<5 Gyr). 
Among the complete sample of well characterized hot 
Jupiters, HATS-9b and CoRoT-17 b (10.7 ± 1.0 Gyr) 
are the oldest systems known to have an age uncertainty 
better than 20%. Figure [9] shows the radius as func¬ 
tion of age for hot Jupiters with 0.5-Mj < M v <2Afj 
and orbital period P < 10 days having age uncertain¬ 
ties smaller than 40%. Systems older than 3 Gyr exhibit 
the expected contraction of the envelope through time 
but most of them are systematically more inflated than 
expected from theoretical models of structure and evolu¬ 
tion. By fitting a straight line through the planets with 
ages higher than 3 Gyr we obtain an empirical contrac¬ 
tion function for hot Jupiters: R p = 1.45 — 0.03t, where 
t is the age of the system in Gyr. The difference be¬ 
tween the theoretical function and the empirical relation 
decreases with the age of the system and for the case of 
HATS-9b both functions are consistent with the observed 
values. The proposed empirical relation between the age 
of the system and the radiu s of t he planet shown in Fig¬ 
ure [9| supports the study of Burrows et al. p007) where 
for young giant planets the higher opacity produced by 
heavy elements delays the contraction, while at later ages 
the higher mean molecular weight dominates and leads 
to smaller radii. However, in order to perform a precise 
study of the evolution of the the radii of giant extraso¬ 
lar planets, particular models with the properties of each 
system should be constructed. 

A possible confusing factor in Figures 7, 8 and 9 is 
the assumption of zero albedo and complete heat redis¬ 
tribution. The measurement of secondary transits on 
these systems in different wavelengths will be informa¬ 
tive for explaining the departure of HATS-9b from the 
correlation. A more precise determination of the radius 
of HATS-9b is also required. The somewhat larger uncer¬ 


tainty in the radius is a result of the incomplete photo¬ 
metric follow-up for this system. The errors in the planet 
radius are governed at this point by the light-curve data, 
but future precise measurements of the transit of HATS- 
9b will be able to lower this uncertainty until it becomes 
dominated by the uncertainties on the stellar parameters. 

Future precise RV measurements of HATS-10b are re¬ 
quired to determine a more precise mass of the planet 
and to explain the high jitter measured with FEROS and 
Coralie with respect to Subaru/HDS. One possible expla¬ 
nation may be the presence of another planetary com¬ 
panion. Subaru/HDS observations, which don’t seem to 
show enhanced jitter, were performed in three continu¬ 
ous days, while Coralie and FEROS observations were 
separated by months and in this case the influence of a 
second more distant companion should be stronger. The 
jitter values quoted in Table [6] refer to RV uncertainties 
for each instrument that have to be added in quadrature 
to the formal RV errors in order for them to be consistent 
with the RV signal computed with the orbital parameters 
of the system. 


4.1. K2 possibilities 

Even thought HATS-9b and HATS-lOb are located in 
the nominal coordinates of field 7 of K2, only HATS-9b 
falls on working silicon. A proposal to observe this star 
in short cadence was recently submitted. The high pho¬ 
tometric precision of K2 will allow us to estimate a much 
more precise radius for HATS-9b, which will help us in 
determining if this planet is a true outlier in the correla¬ 
tion between planet radius, equilibrium temperature and 
planet mass. The high insolation of this planet makes it 
a very good target for measuring secondary transits and 
phase curve variations with K2, which will allow us to 
estimate the albedo and provide a more reliable estimate 
of its equilibrium temperature. Figure [To] shows a mea¬ 
sure of the reflected light signature , (Rp/a ) 2 , for hot 
Jupiters observed by Kepler as a function of planetary 
radius. From this Figure we can see that the potential 
of detecting reflected light signatures of HATS-9b is high 
and its amplitude should be similar to the one of the 
giant planets observed by Kepler so far. Other subtle 
photometric effects, like ellipsoidal variations, Doppler 
beaming and the measurement of asteroseismological fre¬ 
quencies, if present, will also be very valuable for the 
detailed characterization of this particular planet. 
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Figure 7. (Left): Mass-radius diagram of giant TEPs. HATS-9b is marked with a filled square and HATS-lOb with a fille d trian gle. 
Isodensity curves are plotted with dashed lines for pp = (0.25, 0.5,1.0, 2.0, 4.0} gr cm 3 and the 4.5 Gyr isochrones (Fortney et al.|2007) for 
core masses of 0 and 100 Mj with solid lines. (Right): Equilibrium temperature versus radius diagram for giant TEPs. Again, HATS-9b 
is marked with a filled square and HATS-lOb with a filled triangle. 
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Figure 8. Mass-radius diagram of giant TEPs having similar insolation levels to HATS-9b (1750K< T eg <1900K). HATS-9b is marked 
with a triangle. Filled symbols are colored according to the metallicity of the host star. HATS-9b does not follow the correlation formed 
by the other hot Jupiters with similar irradiation levels. 
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Figure 9. Radius as function of the age of the system for hot Jupiters having 0.5 M,t < M n < 2Mj, P<10 days and an age estimation with 
a precision better than 40%. Green lines are the theoretical models of|Fortney et al.|(|2007|) for M p = lMj, a = 0.02AU and a core mass 
of 0 (dashed) and 50 (solid) times the mass of the earth. The red line is an empirical relation computed with these data points. HATS-9b 
is marked with a triangle. Hot Jupiters older than 3 Gyr follow the contraction of their radius over time but the observed contraction rate 
is steeper than the one predicted from the theoretical models. The theoretical radii for Hot Jupiters with ages greater than 10 Gyr (like 
HATS-9b) is consistent with the observations. 
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Figure 10. Reflected light signature as function of the planetary radius for the hot Jupiters observed with Kepler. The symbols are 
colored according to the planetary mass. HATS-9b is marked with a triangle. Given that the photometric precision of K2 is similar to the 
one of the original Kepler mission, phase curve variations and secondary the secondary transit of HATS-9b should be measured by K2. 
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